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I. INTRODUCTION 

In the conventional approach, the inclusive J/ip productions in various high energy pro- 
cesses have been studied in the framework of perturbative QCD (PQCD) and the color 
singlet model [|l|. In this model, one assumes that the color-singlet cc state is produced in 
the ^5*1 state in the short distance scale {i.e. via some hard process which is calculable in 
PQCD) from the beginning, and this object eventually evolves into the physical J /ip state in 
the long distance scale. However, this approach, when applied to the J/ip or ip production 
at the Tevatron, severely underestimates the productions rate 0. In order to reconcile the 
data and PQCD predictions, a new mechanism for heavy quarkonium productions has been 
suggested 0, the color-octet gluon fragmentation into J/ip. Also, the color-octet mech- 
anism in heavy quarkonium productions at hadron colliders through the color-octet [ccjs 
pair in various partial wave states "^^^^Lj has been considered beyond the color-octet gluon 
fragmentation approach 0], In Refs. [Q and [Q, a large class of color-octet diagrams has 
been considered which can contribute to the J/ip production at hadron colliders. Here, the 
basic picture is the foolowing : at the parton level, one can have color-octet cc states with 
various "^^^^Lj, 

qq^(cc)CS?), (1.1) 
gg^icc)Cs!,'^ or ^pP), (1.2) 

at the short distance scale, and the subsequent evolution of the (cc)8(^'^+^Lj) object into 
a physical J/ip by absorbing/emitting soft gluons at the long distance scale [| The short 
distance process can be calculated using PQCD in powers of a^, whereas the long distance 
part is treated as a new parameter (O|0g (^'^■^^Lj)|0) which characterizes the probability 
that the color-octet {cc)(^^~^^Lj) state evolves into a physical J/ip by emitting/absorbing 
soft gluons 0. By fitting the J/ip production at the Tevatron using the usual color-singlet 
production and the cascades from Xc(l-P) and the color-octet contribution, the authors of 
Ref. determined 

{0\Ot{^Si)\0) = (6.6 ± 2.1) X 10-3 GeV^ 

mipm + = (2.2 ± 0.5) . 10- Gev^ 

for Mc = 1.48 GeV. Although the numerical values of two matrix elements {0\Og{^Po)\0) 
and (0 1 (^5*0)10) are not separately known in Eq. (1.4), one can still extract some useful 
information from it. Since both of the color octet matrix elements in Eq. (1.4) are positive 
definite, one has 



^ This kind of processes is clearly possible, but has never been addressed before in a quantitative 
way. The color evaporation model Q is in principle close to this picture, but there are some 
parameters which are introduced phenomenologically, and thus have no definitions in QCD. 

^See Ref. for more details on NRQCD, and definitions of Og . 
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(1.3) 
(1.4) 



< (0|C»^(^^o)|0) < (6.6 ± 1.5) X 10-2 GeV^ (1.5) 
< ^°'^^^^f°^'°^ < (2.2 ± 0.5) X 10-2 GeV^ (1.6) 

These inequalities can provide us with some predictions on various quantities related with 
inclusive J/ip productions in other high energy processes, which enables us to test the idea 
of color-octet mechanism. 

Since the color-octet mechanism in heavy quarkonium production is a new idea proposed 
in order to resolve the ijj anomaly at the Tevatron, it is important to test this idea in 
other high energy processes with inclusive heavy quarkonium productions. Up to now, the 
following processes have been considered : J/ ip production at the Tevatron and fixed target 
experiments Q p, spin alignment of the color-octet produced J/^/^ 0, the polar angle 
distribution of the J in the e^e~ annihilations into J/ip + X [|l^], inclusive J /ip production 
in B meson decays [0], the Z° decays at LEP [jl3|, and J/ip photoproductions at the 



fixed target experiments as well as at HERA |]T4[ |jT5[ |T6|. These processes also depend on 



the aforementioned three color-octet matrix elements in different combinations from (1.4). 
Thus, one can check if the color-octet mechanism provides reasonable agreements between 
PQCD and the experimental data on inclusive J/ip production rates from these processes. 

In this talk, I review two applications of the idea of the color-octet J/ip productions, 
which I was working on, among many recent applications mentioned above 0. First in 
Sec. |I|, I discuss the B ^ J /ip + X using the factorization formula derived in Ref . |Tl| , and 
find that the relations (1.5) and (1.6) overestimate the branching ratio for B —>■ J/ip + X , 
especially for (1.6). Then, we discuss the J/tp photoproductions in the in Sec. |ITT[ For this 
process, the singlet contribution from j + g —>■ J/iIj + g (the 'jg fusion) has long been known. 
And we consider the color-octet subprocesses 

7 + ^7-(cc)(i5f or3pt)o,2), (1-7) 

which have not been included in previous studies. These color-octet 2—^1 subprocesses can 
also contribute to the 2^2 subprocesses through 

j + g-.(cc)Csj,'^ or 'PP)+g, (1.8) 
7 + icc)CS^'^ or' PP) + q. (1.9) 
^ + g^(cc)CS+'Lf)+g. (1.10) 



All of these color-octet 2 — 2 subprocesses are calculated in Refs. |T^. Numerical 
analyses relevant to the fixed target experiments and HERA have been performed. We show 
that the relations (1.5) and (1.6) yields too large a cross section for the J/ip photoproduc- 
tion in the forward direction. They also leads to too rapidly growing da/dz distribution 
for high z region compared to the experimental observations. All of these seem to indicate 
that the relations (1.3) and (1.4), especially the latter, are probably overestimated by an 



^ The case of J/ip productions at the Tevatron is covered in detail by Cho and Leibovich in 
Refs. [4,5]. 



3 



order of magnitude. This is not surprising at all, since the analyses in Ref. employed 
the leading order calculations for the color-singlet parton subprocess for the J / ip hadropro- 
duction. We summarize our review and speculate the origins of these overestimates of J/i/j 
photoproductions and B meson decays in Sec. 

II. INCLUSIVE J/xp PRODUCTION IN B DECAYS 



The effective Hamiltonian for b — > ccq ( with q = d,s) is written as |T7 

■2C, - C_ 



Gf 

Heff = VcbV*g 



C7m(1 -75)c qY{^-l5)b 



+ (C+ + C7^(l - 75)T"c q^il - 75)T"6] , (2.1) 

where C-t's are the Wilson coefficients at the scale fi ~ M^. We have neglected penguin 
operators, since their Wilson coefficients are small and thus they are irrelevant to our case. 
To leading order in ^^(Mb) and to all orders in as{Mi,) hi{MiY / M^) , the above Wilson 
coefficients are 

C+{Mh) ^ 0.87, C4Mh) ^ 1.34. (2.2) 
According to the factorization theorem for the S*— wave charmonium productions in B 



decays, one has [|T7 



Tib W + X)= (OlQi"^"^^^)!") f ,(6 (cc),(^50 + X), (2.3) 

in the nonrelativistic limit, where Fi are rates for hard subprocesses of b quark decaying 
into a cc pair with suitable angular momentum and vanishing relative momentum in the 
color-singlet : 

Ub ^ icc)^CS^) + s,d) = (2C+ - C^f (^1 + fo, (2.4) 

fi(6 ^ (ccU'So) + s,d) = (2C7+ - C_)2 fo, (2.5) 



with 



r„ . IV-.r f ^) (l - ^ 1 . (2.6) 



The operator Oi (^^'^^Sj) is defined in terms of heavy quark field operators in NRQCD 
^ Its matrix element {O\0i {^^^^Sj)\0) contains the nonperturbative effects in the heavy 
quarkonium production processes, and is proportional to the probability that a cc in a 



^We follow the notations in Ref. and will not give explicit forms for these dimension-six 
operators in this paper. 



4 



color-singlet 5*— wave state fragments into a color-singlet S*— wave cc bound state such as 
a physical J/ip., or ip . It is also related to the matrix element {H\Oi{^^~^^Sj)\H) and the 
nonrelativistic quarkonium wavefunction as follows : 

{0\Oi/^CS,)\0) ^ 3 {JMO,eS,)\J/ij) ^ (|-) \R40)f , (2.7) 

in the nonrelativistic limit. Note that dependence on the radial quantum numbers n enters 
through the nonperturbative parameters, (0|Of^(^S'i)|0). 

Using the leptonic decay width of J/ip and ip , one can determine 

{J/^\OiCSi)\J/tlj) ^ 2.4 X 10"^ GeV^ (2.8) 
{^'\OiCSi)\tlj') ^ 9.7 X 10-2 GeV^ (2.9) 

in the nonrelativistic limit with as{Mc) = 0.27. ^ From these expressions with Mh ^ 5.3 
GeV, one can estimate the branching ratios for B decays into J/tp + X and ip + X : 

5(5 ^ J/^ + X) = 0.23%, (0.80 ±0.08)%, (2.10) 
B{B ^ tp' + X) = 0.08%. (0.34 ±0.04 ±0.03)%. (2.11) 

The recent data from CLEO [l^ are shown in the parentheses, where the cascades from 



B — >• Xcj(l-P) ± X followed by Xcj J/ip + 'y have been subtracted in the data shown. 
In view of these results, we may conclude there are some important pieces missing in the 
calculations of decay rates for B — > (cc)i('^S'i) ± X using the color-singlet model in the 
nonrelativistic limit. 

In view of this, we first estimate the color-octet contributions to — *• J/ip+X, motivated 
by the suggestion that the color-octet mechanism might be the solution to the ip puzzle at 
the Tevatron. Although it is of higher order in f ^ (~ 0(f ^)), it can be important in the case 
of the inclusive B decays into J/ip + X, since the Wilson coefficient of the color-singlet part 
is suppressed compared to that of the color-octet part by a factor of ~ as- (In Eq. (2.1), 
(2C+ — C_) ^ 0.4, and (C+ ± C_) ^ 2.20. ) In Ref. [0, a new factorization formula is 
derived for B ^ J/ip + X : 

v(B TN,-,x) - (o|p//'^(3go|o) \ r )^(^^ ^^^n f 

TiB^J/i; + X)-\ j (2C7+-C_) + ^ ) To 

)J/i'{3Q\\n\ / o/i^2\ 



mri'sm ,^ , ^ , 8m| 



^The radiative corrections in has not been included here for consistency. To be consistent with 
the velocity counting rules in the NRQCD in the Coulomb gauge for the heavy quarkonia [0, one 
has to include the relativistic corrections as well, since v ~ as[Mv) in heavy quarkonium system. 
If one includes the 0{a^ radiative corrections to J/V' — > without relativistic corrections, one 
gets a larger (0|O/'''^(35i)|0) compared to the lowest order result, Eq. (2.8) : (0|O/'''^(3S'i)|0) k. 
4.14 X 10^^ GeV^. Relativistic corrections gives a further enhancement. See Ref. [11] for further 
details. 
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^3(0KfeM,c,,C_y^f. ,2.12, 

with To defined in (2.6). Using the relations (1.3) and (1.4), we estimate the above branching 
ratio to be (for as(M|) = 0.28 in Ref. ^) 

(0.42% X 12.8) < B{B J/ij + X) < (0.42% x 13.8) (2.13) 

which is larger than the recent CLEO data by an order of magnitude 0: 

Be^p{B J/ij + X) = (0.80 ± 0.08)%. (2.14) 

Here, the factor 0.42% in Eq. (2.13) comes from the color-singlet and the color-octet ^5"! 
contributions. Other factor comes from the color-octet ^Sq and ^Pj states, which are very 
large if one assumes the relations (1.5) and (1.6). The situation is the same for B ^ ip +X. 
This is problematic, unless this large color-octet contributions are canceled by the color- 



singlet contributions of higher order in 0{as) which were not included in Ref. |T^. If 
there are no such fortuitous cancelations among various color-octet and the color-singlet 
contributions, this disaster could be attributed to the relation (1.4) being too large compared 
to the naive velocity scaling rule in NRQCD, as noticed in Ref. P]. It seems to be crucial 
to include the higher order corrections of 0(af) for the color-singlet J/ip productions at the 
Tevatron, which is still lacking in the literature. 



III. J/ij PHOTOPRODUCTION 



The inelastic J/?/^— photoproduction has long been studied in the framework of PQCD 
and the color-singlet model [jl9[ pO|. The lowest order subprocess at the parton level for 
'j+p—>-J/ip + X is the 7— gluon fusion at the short distance scale 

j + g^(cc)(^S['^)+g, (3.1) 
followed by the long distance process 

(cc)(35f)) ^ j/^, (3.2) 

at the order of O^aa^v^) in the nonrelativistic limit. Thus, the production cross sec- 
tion is proportional to the gluon distribution inside the proton. This is why the 



^Even if we use the new determination (3.38) by Fleming et al. |15], we still get a large branching 
ratio : 

(0.42% X 3.45) < B{B J/^j + X) < (0.42% x 5.45), 
although the discrepancy gets milder than the case (2.13). 
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J/?/'— photoproduction has been advocated as a clean probe for the gluon structure function 
of a proton in the color-singlet model. Without further details, we show the lowest order 
color-singlet contribution to J/ip photoproduction through 7— gluon fusion in the nonrela- 
tivistic limit : 



AM] 



2\2 



t^{t - m^Y + u{u - AM] 



2\2 



{§ - 4M2)2(t - 4M2)2(u - 4M2)2 



where 



lab 



Z = 

s = [k + q^f 
i = (P - kf 

The overall normalization Mi is defined as 



PN-k ' 



S2 

y 

The parameter Gi{J/ip), which is defined as 



M2 



(3.3) 



(3.4) 



(3.5) 



(3.6) 



3c(l) 



in the NRQCD, is proportional to the probability that a color-singlet cc pair in the -'SI 
state to form a physical J/ip state. It is related with the leptonic decay via 



r(J/V^ ^ l+r) = - T^elo? GiiJ/il)), 

o 



(3.7) 



to the lowest order in as- From the measured leptonic decay rate of J/ip., one can extract 

Gi{J/%lj) ^ 106 MeV, (3.8) 

Including the radiative corrections of 0{as) with as{Mc) = 0.27, it is increased to ^ 184 
MeV. Relativistic corrections tend to increase GilJ/tp) further to ~ 195 MeV ||TT |. 
The partonic cross section for 'j + a—^J/ip + b is given by 



di 167rs2 
The double differential cross section is 



dzdP^ 



{7 + p{pN)^J/i^{P,e)+X) 



xg{x,Q^ 



z(l — z) Wits'" 



where 



s 

X = - 
s 



1 

zs 



(3.9) 



(3.10) 



(3.11) 
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One has the following constraints for x, z, t and P^f : 



Ml 

-^<x<l, (3.12) 
s 

-(S - Mj) < £(= i) < 0, (3.13) 



The z and distributions can be obtained in the following manner : 



dz Jo dzdP^ ^' 

da /-^max d'^cr 



dP^, (3.15) 



dz, (3.16) 



W = ^is + Ml + J{s-Mlf-AsP^) , (3.17) 



2s 



^(^ + ^^-\/(^-^|)'-4sP|). (3.18) 



There are two kinds of corrections to the lowest order result in the color-singlet model 
(3.3) : the relativistic corrections of 0{v'^) and the PQCD radiative corrections of 0{as) 
relative to the lowest order result shown in (3.3). We briefly summarize both types of 
corrections here, since they have to be included in principle for consistency, when one includes 
the color-octet mechanism in many cases. 

The relativistic corrections to the 7— gluon fusion were studied by Jung et al. PU 



They found that relativistic corrections are important for high z > 0.9 at EMC energy 
{y/s^p ~ 14.7 GeV). Since it mainly affects the high z region only, we neglect the relativistic 
corrections, keeping in mind that it enhances the cross section at large z > 0.9. 

The radiative corrections to the J/ip photoproduction is rather important in practice. 
This calculation has been done recently in Ref. [^], and the scale dependence of the lowest 
order result {Q^ in the structure function in Eq. (3.10) ) becomes considerably reduced. 
For EMC energy region, the K factor is rather large, K 2. For HERA, it depends on 
the cuts in z and P^. We include the radiative corrections in terms of a factor suitable 
to the energy range we consider. Another consequence of the radiative corrections to the 
color-singlet J/ip photoproduction is that the PQCD becomes out of control for z > 0.9 
at EMC energy. For HERA, one gets reasonable results in PQCD when one imposes the 
following cuts in z and P^ : z < 0.8 and P^ > 1 GeV^ . Thus, it does not make much sense 
to talk about the z or px distributions for such z region in PQCD. One has to introduce 
cuts in z as well as in pt- Following the Ref. PT|, we adopt the following sets of cuts : 



z<0.9, for EMC, (3.19) 
0.2<z < 0.8 for HERA. (3.20) 

At HERA energies, the lower cut in z{z > 0.2) is employed in order to reduce backgrounds 
from the resolved photon process and the b decays into J/ip. For these cuts, the K factor is 
approximately ~ 1.8 both at HERA and the fixed target experiments. We include these 
radiative corrections to the subprocess (3.1) by setting K ~ 1.8. 
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Let us consider color-octet contributions to the 2^1 subprocesses via 



(3.21) 



followed by {cc)s fragmenting into J/ip with emission of soft gluons. This subprocess occurs 
at O(aasf^). Here, a,b are color indices for the initial gluon and the final color-octet cc 
state, and we are interested mS = L = J = and S = L = 1, J = 0,1,2. There are 
2 diagrams representing the vertex. Here we consider the process where only the gluon is 
off-shell. Following the conventions adopted in the previous section, we first write the matrix 
O related to this effective vertex. 



0{P,q) 



V2 



ab 



7^ 



;f + g - A;)2 - M2 



a 2 



+ 4- ^ + M, 



:i + q-gy-M^ 



(3.22) 



With this matrix O we can derive the effective vertices for the 'yg{cc 



7^N2S+l r (8) 



as 



A<'(35f )) = 0, 



M'CPi""^) 



266,(7^5"" g 



l2Ml 



v^eecgsS'' 



g"^ + 2 



g^ - 4M2 



(3.23) 
(3.24) 
(3.25) 



M2(c/2 _4M2) 



C/2 - 4M2 



ea{J.)elelPr 



M'eP^''>] 



'^-^^^''\m. (g^^g^^ + ,,. ^v- + p-.;;-^v^ \ ^^^^j^^^^^ 



{g'-my 



^2 - 4M2 



(3.26) 
(3.27) 



Since J/iIj can be produced via the 2^1 subprocesses with these effective vertices, we can 
obtain the 2-^1 color octet contribution by using the following average squared amplitudes 
asQ 

)r = 2(eee^7. 
' " = 0, 



is) 1 



6 

Ml 
0, 



M2 



[eccgs 



[eccgs 



(3.28) 
(3.29) 

(3.30) 

(3.31) 
(3.32) 



^Our results agree with those obtained in Refs. [Q 
the invariant matrix is different from theirs. 



. Note, however, that our convention of 



9 



The J /if) photoproduction cross section via 2^1 process is given by 

- 64M5 P^VpWJ,=4iV4V. M2 + 7 J' ^^'^"^^ 

Since o" oc 5(1 — z), this 2^1 color-octet subprocesses contribute to the elastic peak in the 
J/'?/'— photoproduction. It is timely to recall that the color-singlet model with relativistic 



corrections still underestimates the cross section for 2; > 0.9 by an appreciable amount pO 
As z — >• 1, the final state gluon in the 7— gluon fusion becomes softer and softer, although 
this does not cause any infrared divergence in the transition matrix element. Therefore, 
it would be more meaningful to factorize the effect of this final soft gluon into the color- 
octet matrix elements, (Og (^^o)) and {Og(^Pj)). The color-octet ^5*0 and ^Pj states might 
reduce the gap between the color-singlet prediction and the experimental value of da/dz for 
0.9<z< 1. 

The color-octet 2^1 subprocess (3.21) considered before not only contributes to the 
elastic peak of the J/ip photoproduction, but it also contributes to the resolved photon 
processes at O(aaff^), where the initial partons can be either gluon or light quarks [q = 
u, d, s). These processes are suppressed by v"^ but enhanced by l/a^, relative to the resolved 
photon process in the color-singlet model. Also, they can enhance the high-pT^J/^/^'s, which 
might be relevant to the J/ip photoproduction at HERA. This can be a background to the 
determination of gluon distribution function of a proton, if the cross section is appreciable. 
The resolved photon process in the color-singlet model is dominant over the 7— gluon fusion 
in the lower z region, z < 0.2, and it can be discarded by a suitable cut on z. Since the color- 
octet contribution to the resolved photon process has not been studied in the literature, we 
address this issue here. When one considers 

j + g^{cc)sCSoOT 'Pj), 
j + q^icc)sCSoOT 'Pj), 

one has to include 

^ + g-^{cc)sCS,) 

simultaneously, since both are the same order of O^aa^v"^). This diagram is the same as the 
color-singlet case except for the color factor of the (cc) state. 

It is straightforward, although lengthy, to calculate the amplitudes for the above three 
processes. Using REDUCE in order to the spinor algebra in a symbolic manner, we can get 
the averaged squared for various 2 — >• 2 processes. 

Another color-octet {cc){^Si'^) contribution to the J/V'— photoproduction comes from 
the Compton scattering type subprocesses : 

l{k,e)+qip,) {cc)est^)iP,en + q{p2), (3.34) 

where P and e* are the four momentum and the polarization vector of the ^Si color-octet 
state, and a is its color index. This subprocess, if important, can be a background to the 
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determination of the gluon distribution function in a proton, since it is initiated by light 
quarks. From the naive power counting, however, we infer this subprocess occurs at O(aa^) 
in the couphng constant expansion, and also suppressed by compared to the color-singlet 
contribution (3.1) due to its color-octet nature. Thus, this subprocess is expected to be 
negligible. 

One can actually quantify this argument by explicitly evaluating the FejTiman diagrams 
for (3.34). The effective vertex for qq {cc){^S[ ) is given by |^ 

M\qipMP2) - (cc)('5r)) = ^ vip2)rT'^uip,) e;ip^+p2,S.), (3.35) 

where e* is the polarization of the produced spin-1 color octet object. Using this effective 
vertex, one can calculate the amplitude for 75 {cc)s{^Si^^"'^)q 



M\jq^icc)CSr)q) = -^-^u{p,) 



u{pi). (3.36) 



where eCg is the electric charge of the light quark inside proton(g = u,d,s). The average 
amplitude squared for the color-octet ^5*1 state is given by 

ElM'iiq - {C-C)CS[%)\' = -:^i9'see,rd + | + 8^). (3.37) 

This completes our discussions on the color-octet 2^2 subprocess for J/ip photoproduc- 
tions. 

Now, we are ready to show the numerical results using the analytic expressions obtained 
in the previous section. Let us first summarize the input parameters and the structure 
functions we will use in the following. The results are quite sensitive to the numerical values 
of Us and and the factorization scale Q. We shall use as(M^) = 0.3, = 1.48 GeV 
and = (2mc)^. For the structure functions, we use the most recent ones, MRS A p2| and 



CTEQ3M , which incorporate the new data from HERA , on the lepton asymmetry 
in boson production [^] and on the difference in Drell-Yan cross sections from proton 
and neutron targets p6[. For the 2 — > 1, we show results using both structure functions. 
For the 2-^2 case, we show the results with the CTEQ3M structure functions only, since 
the MRSA structure functions yield almost the same results within ~ 10% or so. 

Let us first consider the J/ip photoproduction via the color-octet 2 — 1 subprocess. 
Since the subprocess cross section (3.33) vanishes except at 2; = 1, one can infer that it 
contributes to the J/ip photoproductions in the forward direction {z ~ 1,P^ — 0). In 
Figs. 1 (a) and (b), we show the J/ip photoproduction cross section in the forward direction 
{(Jforward) as wcll as the data from the fixed target experiments and the preliminary data 
from HI at HERA, respectively. In each case, the upper and the lower curves define the 
region allowed by the relation (1.4) for two color-octet matrix elements, (O|0g (^5*1)10) and 
(OjOg (^Po)|0). In case of fixed target experiments, (Jforward is usually characterized by z > 
0.9, with the remainder being associated with the inelastic J/ijj photoproduction. According 
to this criterion, the experimental value of 0"exp(7 + p —>■ J/ip + X) contains contributions 
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from inelastic production of J/ip^s. Thus, the data should lie above the predictions from 
the color-octet 2 — >• 1 subprocess, (3.21). Fig. 1 (a) shows that the situation is opposite to 
this expectation. Color-octet contributions are larger than the data, which indicates that 
the numerical values of the color-octet matrix elements are probably too large. At HERA, 
one has the elastic J/ip photoproduction data, which can be identified with the color-octet 
2 1 subprocess. By saturating the relation (1.4) by either color-octet matrix element, we 
get the J/ip photoproduction cross section in the forward direction (Fig. 1 (b)). We observe 
again that the color-octet contribution with (1.4) overestimates the cross section by a large 
amount. This disagreement can arise from two sources : (i) the radiative corrections to 
pp J/tp + X, which were ignored in Ref. is important, and/or (ii) the heavy quark spin 
symmetry for (O|0g (^Pj)|0) ~ (2 J + 1) (OlOg (^Po)|0) may not be a good approximation. 
Although the heavy quark spin symmetry relation is used quite often in heavy quarkonium 
physics, it may be violated by a considerable amount fl^ . 



Recently, Amundson et al. performed the fit to the available fixed target experiments 
and the HERA data independently, and found that [13] 

{Q\Ot{^S,m + ^ (0|08^(3Po)|0) = (0.020 ± 0.001) GeV^ (3.38) 

using the MRSA^'), and CTEQ3M structure functions with as{2M^) = 0.26 and = 
1.5 GeV. This determination is not compatible with the relation (1.4), since the resulting 
(Osi^Po)) is negative. This is another way to say that the determination of the color-octet 
matrix elements from the J/ip productions at the Tevatron may not be that reliable. In fact, 
this is not very surprising, since the radiative corrections to the lowest-order color-singlet 
contributions to the J/ip hadroproductions are not included yet. 

Next, we consider the J/ip photoproduction through 2 — 2 parton-level subprocesses. 
As discussed before, the PQCD corrections to the lowest order ■y + g J/ip + g is not under 
proper control for z > 0.9. Therefore, we impose a cut z < 0.9 at EMC energy, ^ys^ = 14.7 
GeV, and at HERA with ^Js^p = 100 GeV, we impose cuts on z and ||2l]] : 



0.2< z < 0.8, P| > 1 GeV^. 

In both cases, we set K ~ 1.8. 

In Figs. 2 (a) and (b), we show the da/dz distributions of J/ip at EMC (NMC) and 
HERA along with the corresponding data. In both cases, the color-octet ^Si contribution 
(Compton scattering type) is negligible in most regions of z, and thus can be safely neglected. 
The thick dashed and the thin dashed curves correspond to the cases where the relation (1.4) 
is saturated by (0 1 Og ('^Pj) |0) and (O|0g (^S'o)|0), respectively. The thick and the thin solid 
curves represent the sum of the color-singlet and the color-octet contributions, in case that 
the relation (1.4) is saturated by {0\OtCSo)\0) and {0\Ot Po)\0) , respectively. In either 
case, we observe that the color-octet ^5*0 and ^Pj contributions begin to dominate the color- 
singlet contributions for z > 0.6, and become too large for high z region considering we 
have not added the enhancements at high z due to the relativistic corrections. Thus, this 
behavior of rapid growing at high z does not agree with the data points at EMC and HERA, 
if we adopt the determination (1.4) by Cho and Leibovich [^. 

In Fig. 3, we show the inelastic J/ip photoproduction cross section as a function of ^/s^ 
with the cut, z < 0.8 and P^ > 1 GeV^. Again, the color-octet ^5*1 contribution is too small. 
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and thus not shown in the figure. Here again, the color-octet ^5*0 and ^Pj contributions via 
2 — s> 2 subprocesses dominate the color-singlet contribution, if the relation (1.4) is imposed. 
Although the total seems to be in reasonable agreement with the preliminary HI data, 
direct comparison may be meaningful only if the cascade J/'^'s from h decays have been 
subtracted out. There are also considerable amount of uncertainties coming from Mc and 
ag. Therefore, it is sufficient to say that the color-octet ^5*0 and ^Pj states dominate the 
singlet contribution to the photoproduction, if the relation (1.4) is imposed. 

IV. CONCLUSION 

In summary, we considered the color-octet contributions to (A) the inclusive J/ip pro- 
ductions in B decays, and (B) the J/'^ photoproductions (7 + p ^ J/t/j + X) through (i) 
^yg {ccjsi^So and '^Pj) and the subsequent evolution of {ccjg into a physical J/ip with 
z ^ 1 and P^ ^ 0, (ii) the subprocesses 7 + 5'(or q) {cc)g{^So or ^Pj) + g{oT q). These are 
compared with (i) the measured J/ ip photoproduction cross section in the forward direction, 
and (ii) the z distributions of J/ip at EMC and HERA, and the preliminary result on the 
inelastic J/V' photoproduction total cross section at HERA. One finds that the relation (1.4) 
color-octet lead to too large contributions of the color-octet ^5*0 and '^Pj states to the above 
observables. Especially, the first two observables contradict the observation. This is also 
against the naive expectation that the color-octet contribution may not be prominent as in 
the case of the J/ ip hadroproductions, since they are suppressed by (although enhanced 
by one power of a^) relative to the color-singlet contribution. It is also pointed out that the 
same is true of the process B ^ J/ip + X, in which the relation (1.4) predicts its branching 
ratio to be too large by an order of magnitude compared with the data. 

Therefore, one may conclude that the color-octet matrix elements involving ccsC^Sq, ^Pj) 
might be overestimated by an order of magnitude. Since the relation (1.4) has been extracted 
by fitting the J/ t/j production at the Tevatron to the lowest order color-singlet and the color- 
octet contributions, it may be changed when one considers the radiative corrections to the 
lowest order color-singlet contributions. 
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Figure Captions 



Fig. 1(a) The cross sections for '-f + p J/ip + X in the forward direction at the fixed 

target experiments as a function of E^. The solid and the dashed curves were obtained using 

1 (s) 

the CTEQ3M and the MRS A structure functions. Here, TOT^ is the 5*0 saturated curve 
and TOTp is the ^Pj^^ saturated one. 

Fig. 1(b) The cross sections for 'y + p — > J/tp + X in the forward direction at HERA 
as a function of the square root of s^p. The solid and the dashed curves were obtained using 
the CTEQ3M and the MRS A structure functions. Here, TOT^ is the ^Sq'^ saturated curve 
and TOT„ is the ^P?^ saturated one. 
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Fig. 2(a) The differential cross sections da jdz for 7 + p J + X at EMC as a function 
oi z = Eji^/E^. The singlet contributions are in the thick dotted curve, the color-octet ^Sq 
contributions in the thick dashed curve (with (O^(^^o)) = 6.6 X 10-2 GeV^), and the color- 
octet ^Pj contributions in the thin dashed curve (with {Ot{^Pj))/M^ = 2.2 x lO^^ GeV^). 
The total is shown in the solid curve. The relation (1.4) allows the region between two solid 
curves. Here, TOT^ is the ^5*0^^ saturated curve and TOTp is the ^P^f^ saturated one. 

Fig. 2(b) The differential cross sections da/dz for j+p — > J/ip + X at HERA as a function 
of 2; = Ej/^/E^. The singlet contributions are in the thick dotted curve, the color-octet ^Sq 
contributions in the thick dashed curve (with {OtCSo)) = 6.6 X 10-2 GeV=^), and the color- 
octet ^Pj contributions in the thin dashed curve (with {Oti^Pj))/M^ = 2.2 x 10"^ GeV^). 
The total is shown in the solid curve. The relation (1.4) allows the region between two solid 
curves. Here, TOT^ is the ^Sq'^ saturated curve and TOTp is the ^P^ saturated one. 

Fig. 3 Total inelastic J/tjj photoproduction cross section for z < 0.8 as a function of the 
square root of Sjp. The singlet contributions in the thick dotted curve, the color-octet ^Sq 
contributions in the thick dashed curve (with (OtCSo)) = 6.6 x 10"^ GcV^), and the color- 
octet ^Pj contributions in the thin dashed curve (with {Ot{''Pj))/M^ = 2.2 X 10-2 GeV^). 
The total is shown in the solid curve. The relation (1.4) allows the region between two solid 
curves. Here, TOT^ is the ^^q^"* saturated curve and TOTp is the ^Pj^"* saturated one. 
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